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We discuss the net quark and isovector fluctuations as well as off-diagonal quark flavor suscepti- 
bilities along the chiral phase transition line in the Nambu-Jona-Lasinio (NJL) model. The model 
is formulated at non-zero quark and isospin chemical potentials with non- vanishing vector couplings 
in the iso-scalar and iso-vector channels. We study the influence of the quark chemical potential 
on the quark flavor susceptibilities in detail and the dependence of the results on model parameters 
as well as on the quark mass. The NJL model findings are compared with recent lattice results 
' obtained in two-flavor QCD at finite chemical potential. On a qualitative level, the NJL model 

provides a consistent description of the dependence of quark number fluctuations on temperature 
and baryon chemical potential. The phase diagram and the position of the tricritical point in the 
O^l ' NJL model are also discussed for different parameter sets. 
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During recent years the phenomenological importance of fluctuations in finite temper ature/finite density QCD has 
(f-) | been widely recognized [l[ . The analysis of fluctuations is a powerful method for characterizing the thermodynamic 
' properties of a system. In particular, enhanced fluctuations are an essential characteristic of phase transitions 0, [H, 
0, S S 0, 9- Therefore, modifications in the magnitude of fluctuations have been suggested as a phenomenological 
probe of deconfinement and chiral symmetry restoration in heavy- ion collisions [3, Q. In this context, fluctuations 
related to conserved charges, like baryon number, electric charge and isospin, are particularly relevant This is 
due mainly to the difference in mass gap of the effective degrees of freedom in the confined (hadronic) and deconfincd 
(quark-gluon plasma) phases. Thus, a large change in charge fluctuations is to be expected when a system is passing 
' from quark-gluon plasma to hadronic medium . In heavy ion collisions fluctuations may also reveal information 
on expansion dynamics of the medium created in the initial state and its electromagnetic emissivity as the fluctuations 
pi ' of isospin charge are connected with space-like screening limit of the retarded photon self-energy [l(| . 

, A measure of the intrinsic statistical fluctuations in a system close to thermal equilibrium is provided by the 
i-C corresponding susceptibilities \- Hence, fluctuations in a thermodynamic system may be explored, at least in a 
limited scope, by finding the dependence of x on t ne thermal parameters. The properties of a strongly interacting 
statistical system are characterized by the temperature and a set of chemical potentials. The latter are related to the 
conservation laws, which follow from the global symmetries of QCD. For an isospin symmetric and electrically neutral 



system the thermodynamical ensemble depends only on two parameters, the temperature T and the quark chemical 
■ ■ ' potential. 

One of the consequences of QCD is the existence of a phase transition between the confined, chirally broken hadronic 
phase and the deconfined, chirally symmetric quark-gluon plasma phase. The two phases are separated by a phase 
boundary in the (T, ^t g )-plane. The existence of the phase boundary for < (J. q /T < 1 has rece ntly been established 
by first principle Lattice Gauge Theory (LGT) calculations at finite baryon chemical potential [ill . [i~3l , [Tij . 

Arguments based on effective models [ill Ha, [13, HH, Hsl US HH, HH indicate that at large \x q the transition is first 
order. On the other hand, for small /j, q and two massless quark flavors, the chiral transition is expected to be second 
order with the critical exponents of the 0(4) spin model J23|. For finite quark masses, the second order transition 
is, due to explicit chiral symmetry breaking, most probably replaced by a rapid crossover. The different nature of 
the phase transition at low and high \i q suggests that the QCD phase diagram should exhibit a critical end point 
(CEP), where the line of first order phase transitions end. Beyond the critical end point the transition would then be 
continuous, i.e., second order or crossover. The critical properties of the second order critical end point are expected 
to be determined by the Ising model universality class [20j, |24[ . 

The existence of a critical end point in QCD has been recently studied in lattice calculations at non-vanishing 
chemical potential by either considering the location of Lee- Yang zeros in (2+I)-flavor QCD [HI, or by analyzing 



the convergence radius of the Taylor series for the free energy in 2-flavor QCD [l2j,[26|. Recent results [ll| obtained 
within the first approach suggest that a critical endpoint indeed exists and that it is located at T ~ 164 and fi q ~ 120 
MeV. On the other hand, in ref. [2(| no direct evidence for the existence of a critical endpoint has been found in 
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2-flavor QCD with relatively large quark masses, in the range where the Taylor expansion is applicable, i.e., for 
fj, q < T. 

The critical behavior and the position of the critical end point can also be extracted from observables, that reflect 
the singular part of the free energy. Such an observable is the quark susceptibility Xij > defined as the second derivative 
of the thermodynamical potential f2(T, jl, V) with respect to quark-flavor chemical potential [if, 

x "' = -vlw^> (L1) 

where for two light (u,d)-quarks, jl = (fx u , fid,)- In particular, it was recently argued that the net quark number 
susceptibility Xq may be used to identify the chiral critical point in the QCD phase diagram The flavor off- 

diagonal susceptibility was also investigated at finite temperature in perturbative QCD [28[. In two-flavor QCD, for 
an isospin symmetric system, \ q is given by the sum of the uu and ud susceptibilities: \q = 2(x«u + Xud)- According 
to universality arguments, Xq diverges at the critical end point in the chiral limit as well as for finite quark masses, 
while away from the critical end point, the net quark number fluctuations are expected to be finite. Consequently, 
one expects a non-monotonic structure in Xq along the phase boundary if the QCD phase diagram features a critical 
end point. 

The quark susceptibilities have been computed also in lattice QCD [l3|, [H, [H, [3(| ■ Recently results for the quark 
number and isovector susceptibilities, or equivalently for the diagonal and off-diagonal quark flavor susceptibilities 
have been obtained on a lattice with two light quark flavors using a p4-improved staggered fermion action with a 
quark mass m q /T = 0.4 [l2l, [26j. Results for finite quark chemical potential were obtained by means of a Taylor series 
expansion. The susceptibilities were calculated up to fourth order in the quark chemical potential [26l . l3ll . |32| . 

In general the net quark number susceptibility Xq shows a strong suppression of the corresponding fluctuations in 
the confined phase and a strong increase with temperature near the transition temperature To. Furthermore, the 
susceptibility near the transition temperature To shows a strong increase with increasing quark chemical potential, 
leading to cusp-like structure at \x q ~ T. Besides, the lattice results confirm the expectation that the isovector 
quark susceptibility xi does not exhibit a peaked structure near the transition temperature and shows rather weak 
dependence on \i q . This behavior of the susceptibilities is consistent with the existence of a singularity of the thermo- 
dynamic potential close to /j, q = T = To, the chiral end point. However, the enhancement of Xq a t finite fi q below the 
transition temperature can also be interpreted in terms of an enhanced contribution from baryon resonances; a very 
good description of the T and [i q dependences of the various quark susceptibilities in the confined phase is provided 
by the resonance gas partition function [H, [3l|, . 

Lattice calculations also show a strong correlation between fluctuations in different flavor components. This is 
particularly clear in the LGT results for the off-diagonal susceptibility Xudj which shows a str ong increase of u—d 
correlations near the transition temperature To and an abrupt loss of correlations just above To [26l. l3ll. |32| . 

In this paper we explore the properties of different quark susceptibilities in terms of an effective chiral model. Of 
particular interest is the characteristics of the quark number susceptibilities in different channels along the phase 
boundary and in the vicinity of the critical end point. The calculations will be done in the two- flavor Nambu-Jona- 
Lasinio (NJL) model [34| formulated at finite temperature and chemical potentials for the baryon number and isospin 
densities. 

The NJL model has been used as a model for exploring qualitative features of the restoration of chiral symmetry 
in QCD HH, IHj H3, Hl|- In the particular case of S\J r (2) color symmetry this model has been recently argued to 
provide a realistic description of some of the lattice results [3{|. However, in the SU C (3) case it is rather unlikely that the 
NJL model can provide a quantitative understanding of LGT thermodynamics, since it does not exhibit confinement. 

Consequently, the hadronic degrees of freedom, in particular the baryonic resonances, which provide a quantitative 
interpretation of the LGT results in the confined phase [3ll. |32|. are missing in the model. Furthermore, there arc no 
gluon degrees of freedom in the NJL model. In QCD thermodynamics, in the chirally symmetric/deconfmed phase, 
the gluons play an essential role. Finally, the model suffers from a strong dependence on the ultra-violet cut-off. 

In the chiral limit, the phase diagram of the NJL model shows a phase separation line, where the spontaneously 
broken chiral symmetry is restored. At small densities the transition is second order, while for an appropriate choice of 
the coupling constants, the transition at large densities is first order. Thus, the model reproduces the gross structure 
of the phase diagram expected for QCD. Consequently, the NJL model, formulated at finite T and jl, can be used to 
explore the qualitative behavior of quark susceptibilities and, more generally, universal features of the chiral phase 
transition in the (T, /x g )-plane. 



^ Recently an interesting extension of the NJL model was proposed that mimics confinement by including the Polyakov line as dynamical 
field which couples to constituent quarks [40ll , 
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The net quark number susceptibility has been computed in the NJL model at finite T some time ago Q and 
recently also at finite quark chemical potential \i q Q ■ Our analysis is going beyond previous studies by extending the 
NJL Lagrangian to finite quark and isospin chemical potential as well as to non-vanishing vector coupling among the 
constituent quarks. This allows us to model the net quark number susceptibility Xqi the isovector one xi or equivalcntly 
the diagonal and off-diagonal quark flavor susceptibilities Xff an d to study their dependence on temperature and 
chemical potential as well as on model parameters. On a qualitative level the results can then be confronted with 
recent lattice results. 

The paper is organized as follows: In section ® we introduce the NJL model Lagrangian and its thermodynamics. 
In section [3j we introduce the flavor diagonal and off-diagonal susceptibilities and calculate their T and /i 9 as well as 
model parameters dependences. We discuss the qualitative comparison of the model results with the recent lattice 
findings. Finally in section [U we give a brief summary and discussion of our results. 

2. THE TWO-FLAVOR NJL MODEL 

For two quark flavors and three colors the Lagrangian of the Nambu-Jona-Lasinio (NJL) model reads (HI l35l. l37l. |38| : 

C = "4){i</) — m)ip + Gs (■tpip) 2 +(ipiTj 5 ip)' 

- (Vhv^) 2 ~ G<7 } [(^r>V) 2 + (V^ 7m 75^) 2 ] + ^7oV> , (2.1) 

where m = diag(m„, to^) is the current quark mass matrix, [i = diag(/i u , fid) the chemical potential matrix and f 
denotes Pauli matrices. The strength of the interaction between the constituent quarks in the scalar and vector chan- 
nels is parameterized by the dimensionful coupling constants Gs, Gy and Gy^. We note that while the strength 
of the scalar-isoscalar and pseudoscalar-isovcctor interactions are equal, due to constraints from chiral symmetry, the 
vector-isoscalar and vector-isovector interaction terms are separately invariant, and hence the corresponding interac- 
tion strengths can be chosen independently. In the following calculations focused on the quark number susceptibility, 
we choose the axial- vector condensate to be zero since it does not couple to the vector current. 

The constraints imposed by the conservation of the net quark number of different flavors are controlled by the 
chemical potential \x = diag(/i u , [id). The chemical potentials for the total net quark density n q and the iso- vector 
quark density nj are obtained as linear combinations of \i u and fid 

= + Hd) , Hi = tj(Mu - Md) • (2.2) 
In terms of fi q and /x/, the last term of the Lagrangian can be expressed as 

Ly, = tpfijoip = fiqtp^ip + Hiip^nip . (2.3) 

In addition to the current quark masses and the three coupling constants introduced above, one additional parameter 
is required to complete the model. This is the momentum cut-off (A), which regulates the ultraviolet divergencies. 
In vacuum the values of A and Gs are fixed by requiring that the pion decay constant / ff = 92.4 MeV and the 
pion mass m, = 135 MeV are reproduced. Choosing the current quark masses m u ~ nid = 5 MeV one finds for a 
three-momentum cutoff A = 664.3 MeV and G S A 2 = 2.06 [H|. 

In the chirally broken phase, the ratio of the coupling constants of lo and p mesons to nucleons is empirically 
given by g U NN / 9pNN — 3. This value is also consistent with the naive quark model for the nuclcon, where the 
corresponding couplings to quarks are identical, i.e., g U QQl ' 9pQQ = 1- We account for this on a qualitative level by 
setting Gy^ ~ SG^ in the broken phase and Gy^ = G^ in the symmetric phase. Thus, we consider G^ /Gy^ in 

the range from | to 1, keeping the vector-isoscalar coupling fixed Gy^ — 0.3 Gs- 

The thermodynamics of the NJL model (|2.ip at finite temperature and non vanishing net quark and isospin chemical 
potentials is obtained from the partition function Z(T, \x qi pj, V). In the mean field approximation [l6| the partition 
function is obtained from the effective Lagrangian 

C = ^ -M + ^70)^ - ^tr ((M - m) 2 ) 

+ 4^ (A9_M9)2 + 4G|7 y(A/_/i/)2 ' (2 ' 4) 
where M = diag(M u , Md) and the trace tr is in flavor space. 
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In Eq. (|2.4p we have introduced a dynamical mass M and a shifted chemical potential \x given by 

M = m-2G s {$ip) , 
A = K + ArT3 , 

where 

A« = M<? - 2G V ' (V'70'0) i 
fa = /i/ - 2G|7 } (^370^) • 
The resulting thermodynamic potential density* 2 is of the following form 



(2.5) 
(2.6) 



(2.7) 



w(T,(i;M,p,)= Vf(T,[j,;Mf,(if) 

f=u,d 

1 i 

(A<2 -M<?) 2 



4G S 



tr ((M 



(S 



where 



4G 



u f {T, l i i M f ,iif) = -2N c 



4G 



(V) 



(2.8) 



d 3 p " 



(2tt)3 

,(- 



E/^-TKl-nW^T,^)) 



Tm(l-n ( , \p,T,fa)) 



(2.9) 



In Eq. (|2.9p . Ef(p) = yp 2 + Mj is the scalar part of the quasiparticle energy. The contributions of the vector 

potentials are absorbed in the shifted chemical potential fa. Furthermore, n* is the distribution function for particle 
(+) and anti-particle (— ) states 



(p, T, fa) = ( 1 + exp[(E f (p) T fa)/T] 



(2.10) 

The condensates appearing in Eqs. (|2.6[) - (|2.8[) are determined by extremizing the thermodynamic potential* 3 with 
respect to the dynamical mass and the shifted chemical potentials at a given temperature T and chemical potential p. 



duj duj 



doj 



dAI dft q dfn 



0. 



(2-11) 



These conditions yield the scalar condensate and the quark densities as functions of temperature and chemical po- 
tential. Furthermore, they imply that the scalar and vector fields can be fixed when computing thermodynamic 
derivatives. Thus, one obtains the standard thermodynamic relations for instance for the quark density 



dw(T lf i) 



dfi 



du(T,p;M f ,fa) 



dfi 



T;M,ft 



(2.12) 



However, as discussed below, the dependence of AI and \i on temperature and chemical potential yields non-trivial 
contributions to second derivatives of the thermodynamic potential, e.g., susceptibilities. 
The stationarity conditions (|2.11[) . together with Eqs. (|2.6[) - (|2.8p . imply 



M f = mf + AGsN c ^ 

f=u,d 



d 3 p Mf 



(2tt)3 E f 



1 - n ( f +) [p, T,fa)- n) ' (p, T, fa) 



p q =p q + 4G { ^N c 

f=u,d 



d 3 p 



n f +) (p, T, fa) - n) ' (f, T, fa) 



(-), 



p I = fa + AG^ r) N c 



d 3 p 
(2^ 



(p, T, fa) - n[ ) (p, T, faj\ - (u -»■ d) 



(2.13) 
(2.14) 
(2.15) 



# 2 The thermodynamic potential is given by Q = luV, where V is the volume of the system. 

* 3 The thermodynamic potential is minimized with respect to variations of the scalar field and maximized with respect to variations of the 
(zeroth component) of the vector fields. 
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FIG. 1: The left-hand figure shows the dynamical u-quark mass M u in the chiral limit as a function of temperature T and 
the quark chemical potential {i q . The right-hand figure represents the quark number density n q in the chiral limit in the 
(T, ^t q )-parameter space. The calculations were done for isospin symmetric matter, i.e. hi = 0, with G v = 0.3 Gs- 



By comparing Eqs. (|2.7p . (|2.14[) and (|2. 15[) . we find explicit expressions for the quark number density n q = (ipjQip) 
and the iso-vector density ri\ = (ipTsjoip) 

n q = 2N C J2 I 7S3 Wf +) & T > M - 4 _) & T ' , (2-16) 



m = 2N C 



f=u,d 

d 3 p 



(2tt) 3 



(2tt) s 

(nl+' (p, T, m«) ~ (p, T, - (« -> d)] ■ (2-17) 



In practice, one first solves the gap equation (|2 . 13[) for fixed T, jl q and p/ and then computes \i q and p/ as well as n q 
and n/ using (|2Ti]) - (|2T7]) . 

In Fig. [T] we show the dynamical quark mass Mf and the net quark number density n q in the (T, p 9 )-plane for 
vanishing iso-vector chemical potential fij = in the chiral limit. For (ij = the isovector density n/ vanishes for 
all values of T and p 9 . In Fig. [2] the phase diagram of the NJL model in the (T, /z g )-plane is shown for an isospin 
symmetric system in the limit of vanishing current quark masses. The boundary between the chirally broken and 
symmetric phases was located by finding the onset of chiral symmetry restoration, M(T,fi q ) = 0, when approaching 
from the broken phase. As discussed in the introduction, the order of chiral phase transition is, in the chiral limit, 
expected to change from second order at low to first order at high net baryon densities. Thus, somewhere along the 
phase boundary one expects a tricritical point (TCP), where the order of the chiral transition changes. Close to the 
phase boundary, the thermodynamic potential, may be expanded in a power series in the order parameter M as in 
Landau theory [4l|: u(T, p 9 ; M) = lo q + \aM 2 + jbM 4 + 0(M 6 ). At a second-order transition a = and b > 0, while 
at a first order one a > and b < 0. (In the latter case, the coefficient of M 6 should be positive for stability.) Thus, 
the tricritical point can be identified by a = b = 0. 

In the NJL model the position of the phase boundary and the TCP depends on the model parameters [38l. l42l, l43l, |44| . 
In Fig. [2] we illustrate the dependence on the vector and scalar couplings Gy and Gs as well as on the momentum 

cut-off A. In the left panel the critical line is shown for different strengths of the vector coupling Gy , keeping G5 

( s) 

and A fixed. With increasing G v , the phase boundary at fixed T is shifted to larger [i q . This is expected, since at 

(S) 

non-zero net baryon density, the vector coupling G v provides a repulsive contribution to the energy of a quark and 
thus to the chemical potential (see Eqs. (|2.7[) . (|2.14p 'h In fact, in the T — jl q plane the position of the line of second 

order phase transitions is independent of Gy^ . This is clear from the fact that the gap equation Eq. (|2.13[) depends 

on G v only through p, q . On the other hand, the location of the first order transition and the position of the TCP 

(S) 

do depend on the vector coupling G v . This dependence is reflected in the shift of the TCP to smaller temperatures 

( s) 

with increasing strength of the vector coupling. We note that for Gy > 0.6 Gs, the transition is everywhere second 
order and there is no TCP. In the right panel of Fig. [2] we show the dependence of the phase boundary on the cut off 
A. We have chosen three parameter sets, summarized in Table HI where the cut off is varied from our standard value, 
A = 664.3 MeV, up to almost 1 GeV, keeping Gy ' — 0. The corresponding values of Gs arc chosen such that all 
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FIG. 2: The NJL model phase diagram in the chiral limit for Gy = 0, 0.3 and 0.6 Gs (left panel) and for the set of parameters 
(a)-(c) of Table. 1 (right panel). The dashed (solid) line shows the location of the second-order (first-order) transition. The 

tricritical point, indicated by a dot (•), is located at (T,fj, g ) = (65,275) MeV for = and at (T,/j, q ) = (42,305) MeV for 



G 



(S) 



0.3 Gs- Both phase diagrams correspond to a vanishing isovector chemical potential, fii = 0. 



A (MeV) G S A 2 

set (a) 664.3 2.060 
set (b) 797.2 1.935 
set (c) 995.5 1.829 



TABLE I: Set of cutoff A and scalar coupling constant Gs used in the model calculations. 



parameter sets yield the same critical temperature at [i„ = 0, i.e. T c = 177 MeV, in agreement with lattice results for 
two-flavor QCD [45|. With increasing A, the phase boundary is again shifted to slightly larger values of the chemical 
potential and the TCP to smaller temperatures. 

In appendix A we give a closed expression for the integrals appearing in the gap equation (|2.13|) at the phase 
boundary. Furthermore, by requiring that the critical temperature at fi = remains fixed, as in the right panel of 
Fig. [21 we obtain a relation between the chemical potential at the T = transition and the cut off A 

*=i&hii 5A - (2 - 18) 

This relation is derived under the assumption that the transition is second order everywhere. Nevertheless, it is well 
satisfied also when the transition is weakly first order at high quark densities. Thus, (|2.18[) provides a quantitative 
interpretation of the shift of the phase boundary found in Fig. [2j 

We conclude that in the limit of vanishing vector coupling, (Gy — * 0), the TCP in the NJL model is located at 

(S) 

T ~ 65 and fj, q ~ 275 MeV. For non-zero G v the position of the TCP moves towards lower temperature and higher 
chemical potential, and for sufficiently large vector coupling, the chiral transition is second-order everywhere. 

The phase diagram in the (T, // 9 )-plane shown in Fig. [2] applies to the chiral limit. A non-zero quark mass in the 
Lagrangian (|2.4p will modify the position of the phase boundary. Furthermore, a finite quark mass breaks the chiral 
symmetry of the Lagrangian explicitly. Consequently, the second order transition at high T and small \i is replaced 
by a cross-over transition and the TCP by a critical end point. The critical behavior of the NJL model for finite 
quark masses is consistent with the results obtained in other effective models [(| [13, [H, H2, IH| ■ 

The position of the phase boundary and the order of the chiral phase transition can be also identified through 
thermodynamic observables. In the Introduction we argued that quark number fluctuations are sensitive probes of 
the phase transition. Furthermore, fluctuations of conserved charges are directly accessible in experiments. Thus, it is 
of interest to explore the behavior of the quark number fluctuations in the vicinity of the phase boundary in effective 
models, like the NJL model. 

In the next sections we formulate quark susceptibilities in the NJL model and explore their dependence on thermal 
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and model parameters. We also consider the influence of finite quark masses on the quark fluctuations and discuss 
the NJL model results in the context of the recent lattice findings. 



3. QUARK NUMBER SUSCEPTIBILITIES 



The net quark number and iso-vector susceptibilities Xq and Xi describe the response of the quark density n q and 
the isovector density ni to the change of the corresponding chemical potentials. Thus, \q and Xi are defined as 
derivatives of n q and ni with respect to \x q and \xi 



Xq 



dn q 
d[i q 



Xi 



dni 

dm 



(3.1) 



The net quark and the isovector densities are in the NJL model given by Eqs. (|2.16[) and (|2.17p . The evaluation of the 
derivatives in (|3.ip . taking the implicit dependence of the dynamical masses Mf and the shifted chemical potentials 
flf on fig and fii into account, yields 



2N C 



f=u,d 



d 3 p 



M f dM 



Ef djlq 



,(-) 



(+h 



Xi 



2N C 
T 



d 3 p 



(27T) 

dfi u 



1 - n 



M„ <)M 



•„W(1 



,(-) 



(3.2) 



«(nW(l-n(+))-n(-)(l-nl-))) 



E u d[ii 

( n (+)(l-4+))+4-)(l-nH))-(t 



d) 



(3.3) 



The 



In Eqs. (|3.2p and p.3p we have suppressed the T and /i/ dependence of the distribution functions 
derivatives of the dynamical masses Mf and the reduced chemical potentials [if entering in Eqs. (|3.2[) and (|3.3[) are 
given in Appendix iBl 

In addition to the fluctuations of the net quark and isovector densities we also introduce the flavor diagonal and 
off-diagonal susceptibilities defined by 



Xff 



Xff 



d 2 LO 



dfifdur 



(3.4) 



with f^f e{u,d}. 

In isospin symmetric matter the susceptibilities Xuu{= Xdd) and \ud are related to Xq and Xi by 



Xv 



^{Xq+Xl), Xud = ^(Xq~Xl) 



(3.5) 



In the following section we compute the susceptibilities introduced in Eqs. ()3.2|) - (|3.5|) in the NJL model and discuss 
the dependence of the quark number fluctuations on temperature and chemical potentials in the vicinity of the phase 
boundary. 



3.1. Quark susceptibilities in the NJL model 



As discussed above, the NJL model does not exhibit the confinement property of QCD. Thus, there are no hadronic 
bound states and resonances in the chirally broken phase. Instead we are dealing with constituent quarks which can 
be viewed as quasi-particles, with a temperature and density dependent mass. At the chiral transition the composition 
of the medium is not changed in the NJL model; the dynamical quark masses Mf vanish and above T c the medium 
is populated by interacting massless quarks. Furthermore, high momentum quark modes are suppressed due to the 
ultra-violet cut-off. This suppression is particularly relevant at high temperatures. The differences in the mass 
spectrum between the NJL model and QCD, as well as the suppression of high-momentum states, results in different 
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FIG. 3: The quark number susceptibility \ q as a function of T/Tq in the chiral limit, for different values of the quark chemical 
potential fi q . Here ' is the quark number susceptibility for an ideal quark gas and To = 177 MeV is the transition 



temperature at ^ q = = 0. The calculations correspond to an isospin symmetric system with G 



(S) 



0.3 G s . 



quantitative properties of the quark number fluctuations. In spite of these differences, the NJL model is useful for 
exploring general features of the susceptibilities close to the phase boundary and, in particular, near the TCP. 

In Fig. [3] we show the quark number susceptibility \ q as a function of T for different values of normalized to 

that of an ideal quark gas Xg frCc) = N c N f (T 2 /3 + fi 2 q /ir 2 ). 

The temperature dependence of \q shows characteristic features, which vary rapidly with \i q . The phase boundary 
is signaled by a discontinuity in the susceptibility. The size of the discontinuity grows with increasing fi q up to 
the TCP, where the susceptibility diverges. Beyond the TCP the discontinuity is again finite. On the other hand, 
at \i q = the discontinuity vanishes and the susceptibility shows a weaker non-analytic structure at the transition 
temperature, corresponding to a discontinuity in dxq/dT. The critical properties of \q are consistent with a second 
order phase transition belonging to the universality class of 0(4) spin model in three dimensions [1, [3l|. Due to the 
lack of confinement in the NJL model, the leading contribution to the thermodynamic potential and to fluctuations 
of conserved charges is due to single-quark loops. The model can be improved by including the interaction of quarks 
with the Polyakov loops. In the resulting PNJL model [U SE|, confinement is mimicked in the sense that three-quark 
states are the leading thermodynamic modes below the phase boundary, in the "confined" phase. In this model the 
fluctuations of the net quark number in the chirally broken phase are suppressed [46| compared to the results of 
the NJL model shown in Fig. 3. Furthermore, in the PNJL model the dependence of the fluctuations on the value 
of the quark chemical potential is also stronger. Recently, it was shown [47| that by an appropriate choice of the 
parameterization of the effective Polyakov loop potential, it is possible to quantitatively reproduce some LGT results 
on fluctuations of conserved charges. Nevertheless, the critical properties of the net quark number fluctuations near 
the phase transition in the NJL and in PNJL models are similar, since both models belong to the same universality 
class. 

We now explore the qualitative features of the critical region within Landau theory [4l|. As already indicated 
in section [2j we construct an effective thermodynamic potential, valid in the vicinity of the chiral transition. The 
thcrmodynamical potential oj(M, T, fi q ) is expanded in a power series in the order parameter M, the dynamical quark 
mass, around M = 0: 

w(T, n q , M) ~ u (T, N ) + i (T, v q )M 2 + h(T, N )M 4 + 0(M 6 ) . (3.6) 

Here we neglect the M 6 term for simplicity. Although it is crucial for the calculation of critical exponents (see 
Appendix [C]) , it does not affect the present argument. We assume that b > 0, i.e., that we are above or at the 
TCP, where the transition is second order. For a > 0, the effective potential (|3.6[) has a minimum at M = 0, which 
corresponds to the symmetric phase, where w(T, /z;0) = u)q(T,h). On the other hand, for a < the minimum is 
located at Mq = y^—a/b and the system is in the broken symmetry phase, where 
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The second-order phase boundary (the 0(4) critical line) is determined by the requirement a = and b > 0. Above 
the critical line, in the symmetric phase, M = and the quark number susceptibility \ q is given by 



The coefficient a(T,fi q ) may be expanded around any point (T c , fj, c ) on the 0(4) critical line. Close to the critical 
line, it is sufficient to keep only the leading terms 

a(T, fi q ) ~ A(T - T c ) + B( Pq - fi c ) , (3.9) 

where the expansion coefficients A and B depend on T c and fi c . 

Using Eqs. (|3.7j) and (|3.9j) obtains the quark susceptibility in the broken phase 

(broken) (sym) , /o iq\ 

Xq Xq + 2&(T,/x g )' K ' 

where we have dropped terms that vanish on the critical line. The discontinuity of Xq across the 0(4) critical 
line at finite \x q is given by the second term in (|3. 10|> . At \x q = the coefficient B, and thus the discontinuity 
of Xqi vanishes by symmetry. Keeping the next term in the expansion of a(T,fj, q ) ~ A(T — T c ) + B2^ q , one finds 
^(broken) ^ = Q ) = Xq^Hftq = 0) + (T — T c )AB 2 /b. Thus, at /x g = the susceptibility at T = T c is continuous, 
while its temperature derivative is discontinuous, as seen in Fig. [3] 

Finally at the TCP both a{T,fi q ) and b(T,/j, q ) vanish. Consequently, the susceptibility in the broken phase (|3.10p 
diverges at the TCP, in agreement with the results shown in the right panel of Fig. [31 With the present choice of 
parameters, the TCP is located at (Ttcp 7 Mtcp) = (42,305) MeV. Beyond the TCP, the phase transition is first 
order. There the susceptibility again exhibits a finite discontinuity at the phase boundary. The susceptibility in the 
symmetric phase, X^ ym \ is expected to vary smoothly along the phase boundary. 

From the perspective of heavy ion experiments, several susceptibilities are of interest. In particular, this applies to 
fluctuations of conserved charges, which may be directly accessible in experiment. Furthermore, those susceptibilities 
that reflect the critical behavior, may possibly be used to explore the QCD phase transition experimentally. As we 
have stressed repeatedly, the quantitative structure of the phase diagram and the position of the critical end point 
are model dependent. Consequently, in detail the QCD phase diagram most likely differs from that found in the 
NJL model. Nevertheless, such a model study can still answer phenomenologically relevant questions concerning e.g. 
the size of the critical region, where the fluctuations are dominated by the singularity at the conjectured critical end 
point. 

In Fig. [5] we show the net quark and isovector susceptibilities \q an d Xi along the phase boundary, given in Fig. [5] 
The position of the TCP is signaled by the singularity of the net quark susceptibility \ q . The corresponding non- 
monotonic behavior of the fluctuations with increasing beam energy may give rise to observable effects in heavy-ion 
collisions [|| # 4 . We find that the critical region, where the fluctuations are dominated by the singularity, corresponds 
to a window AT ~ 30 MeV and Afi ~ 10 MeV around the TCP. In the absence of a TCP, the net quark susceptibility 
would be a monotonic function of T along the phase boundary, as illustrated by the dashed-dotted line in Fig. 2] 
We note that the qualitative behavior of the susceptibility is consistent with the results of Landau theory discussed 
above. First, the discontinuity across the phase boundary vanishes at \x q = 0. Second, the singularity of Xq shows 
up only in the chirally broken phase, while the susceptibility in the symmetric phase is monotonous along the phase 
boundary and shows no singular behavior. 

In nucleus nucleus collisions a change of the collision energy y/s is correlated with a corresponding change of the 
temperature and the chemical potential. An increase of y/s results in an increase of the temperature T and a decrease 
of the baryon chemical potential /i q . Thus, the critical region around tricritical point/critical end point (AT, A[a) can 
be (approximately) converted to a range of center-of-mass energies in A-A collisions. Assuming for simplicity that 
the relation of T c and \fs is the same as for the chemical freezeout parameters extracted from data [48| , we find that 
AT ~ 30 MeV corresponds to Ay/s ~ 1 A-GeV. Consequently, this crude estimate implies that in order to observe 
effects of critical fluctuations in A-A collisions one would need to measure an excitation function with an energy step 
Ay/s smaller than 1 A-GeV. 



* 4 A similar non-monotonic behavior appears in any observable, directly related to the net quark number density-density correlator. Thus, 
measurements of the corresponding non-monotonic structure in the baryon number or electric charge density, net-proton number density 
or in the mean transverse momentum would be excellent experimental probes of the critical end point in the QCD phase diagram. 
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FIG. 4: The quark number (left) and isovector (right) susceptibilities \ q and \i as functions of the temperature along the 
phase boundary. In the left-hand figure the solid (dashed) line denotes \q m the chirally broken (symmetric) phase. The 
vertical dotted-line indicates the position of the tricritical point TCP. The calculations were done in the chiral limit in an 
isospin symmetric system with the vector coupling constant G v = 0.3 Gs- 



In the calculation of the critical properties we have employed the mean-field approximation, which in general does 
not yield the correct critical exponents [49[ . We note however that the non-mean-field critical behavior is suppressed 
near the TCP, since the quartic coupling &(Ttcp, I^tcp) vanishes [f||5(|. Thus, the critical exponents of the TCP are 
close to the mean-field exponents [5l| and mean-field theory provides a good description of e.g. the susceptibilities near 
the TCP. Since the physical quark masses are small, the critical end point is influenced by the tricritical point. Thus, 
one expects the region near the critical end point, where mean-field theory breaks down, to be relatively small [||. 

Both at the TCP and at the critical end point the quark-number susceptibility \ q diverges. However, the critical 
exponents differ. The mean-field exponents of the TCP and the CEP can be obtained from Landau theory. As 
discussed in Appendix [Cl the critical exponent for paths approaching the TCP asymptotically tangential to the phase 
boundary, the susceptibility diverges with the critical exponent 7 9 = 1. Approaching the TCP along the first-order 
transition, the pre-factor is twice as large as along the 0(4) critical line. For other paths the critical exponent is 
lq = \- At the 0(4) critical line, the susceptibility remains finite. The corresponding critical exponent of the 0(4) 
universality class is a ~ —0.2, while in the NJL model we obtain the mean-field value for this critical exponent, a = 0. 
Finally we mention that for non-zero quark mass, at the critical endpoint, the mean-field critical exponent along a 
path not tangential to the phase boundary is 2/3, while along the phase boundary it remains equal to unity Q. When 
fluctuations are included, the first exponent is rcnormalized to that of the 3D Ising model universality class fBll ]. i.e. 
e = 0.78. 

In Fig. [5] we illustrate the critical behavior near the 0(4) critical line and at the TCP. The dependence on the 
reduced temperature t is consistent with the exponents and relative pre-factors obtained in Landau theory. The 
different behavior of the quark number susceptibility at the critical end point and at the 0(4) critical line can be 
traced to the critical behavior of the dynamical quark mass. In Appendix [C] we compute the mean-field exponents 
for the dynamical quark mass in Landau theory. At the 0(4) critical line M 2 ~ \t\, while at the TCP M 2 ~ \t\?. The 
scaling of M 2 obtained in the NJL model are consistent with this, as shown in the right panel of Fig. [H 

Furthermore, in Fig. [6] we also show the "critical" region, where the susceptibility exceeds the free one by more 
than an order of magnitude^ 5 . The different critical exponents along the phase boundary and perpendicular to it 
are reflected in the shape of the critical region. It is elongated along the phase boundary, where the singularity is 
strongest. 

In heavy ion collision an additional complication to explore and map the QCD phase diagram experimentally 
appears due to expansion dynamics, finite system size and secondary hadronic rescattering in a medium. All these 
effects can dilute observation of the critical fluctuations along the chiral phase transition line [52| . 



# 5 By "critical" region wc mean here the region where the susceptibility is large due to fluctuations, not the region of no-mean-ficld critical 
behavior. 
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FIG. 5: The quark number susceptibility near the tricritical point (left) and the scaling of the quark mass as the TCP and 
0(4) critical line is approached (right). At the TCP the reduced temperature is given by t = \T — Ttcp\/Ttcp , while at the 
0(4) line it equals t = \T - T c \/T c . 
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FIG. 6: The "critical" region in the T — n plane obtained in the NJL model. Within this region the susceptibility is enhanced 
by an order of magnitude compared to the free one. 



An interesting observable that characterizes thermal fluctuations related with isospin conservation is the isovector 
susceptibility xi defined in Eq. (|3.ip . The NJL model results for \i m the isospin symmetric system are shown in 
Fig. [7] as a function of T for different The isovector fluctuations, contrary to net quark fluctuations, are neither 
singular nor discontinuous at the chiral phase transition for finite chemical potential. As shown in Fig. |4j we find 
a rather smooth increase of xi with increasing \i q along phase boundary line. At the TCP where the net quark 
number susceptibility diverges, xi remains finite. The non-singular behavior of xi at the TCP is consistent with 
the observation that there is no mixing between isovector excitations and the isosclar sigma field due to SU(2)y 
isospin symmetry [5^]. Also recent LGT results (2(| show a smooth change of the isovector fluctuations around the 
deconfincmcnt transition and a fairly weak dependence of xi on the quark chemical potential [x q . 

The net quark number \q an d the isovector xi susceptibilities are related with fluctuations of the electric charge 
XQ 

1 1 1 d 2 P , 

XQ = m x " + 4 Xl + 6d^- (3 ' n) 
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FIG. 7: The isovector susceptibilities \i f° r different values of net quark chemical potentials /J, q as a function of T/To in the 
chiral limit. The calculations correspond to isospin symmetric system and the vector coupling constant G v = 0.3 Gs- 



Here P is the thermodynamic pressure. For isospin symmetric system the last term in Eq. (|3.11|) vanishes. Hence in 
this case all relevant susceptibilities are linearly dependent. Clearly, since \i is finite at the TCP, the electric charge 
fluctuations XQ diverge with the same critical behavior as \ q . However, at finite [ii the properties of Xi at the chiral 
phase transition in general and at the TCP in particular, change. At non- vanishing fj,j, the SU(2)y symmetry is 
explicitly broken. Thus, the isoscalar sigma field mixes with the isospin density [531 ] - Consequently, the isovector 
susceptibility exhibits a similar structure as Xq-, with a singularity at the TCP. Furthermore, at finite /i/ and away 
from TCP one expects a peak in xi at the chiral phase transition. This is due to the exponential dependence of the 
susceptibility on for massive constituent quarks in the broken phase and the power-law dependence in the chirally 
symmetric phase where the quarks are massless. 



3.2. Model parameter dependence of quark susceptibilities in the chiral limit 

In the previous section, where we considered the net quark and isovector susceptibilities, the coupling constants 
of the effective interaction between constituent quarks were fixed by requiring that the model reproduces vacuum 
observables. In the following we discuss the influence of changes in the model parameters on the critical properties of 
the quark flavor fluctuations. We also present results for the flavor diagonal and off-diagonal susceptibilities Xuu and 
Xud and discuss their properties. 

In Fig. [8] the temperature dependence of the flavor diagonal susceptibility Xuu is shown for several values of the 
quark chemical potential for two choices of the vector couplings Gy . The susceptibility is normalized to the free one 

Xul CC \ which is defined by Xuu C) = Xq 1CC) / N f- 

At vanishing chemical potential there are generic features in the temperature dependence of the quark flavor 
susceptibilities. Consider the net quark Xqi the isovector xi and flavor diagonal Xuu susceptibilities shown in Figs. [3J 
[7]and[5J Clearly all these susceptibilities are strongly enhanced near the chiral phase transition point To, while above 
To the fluctuations are suppressed. The increase of the quark susceptibilities with temperature seen in the broken 
symmetry phase reflects a decrease of the dynamical quark mass as the chiral transition is approached. Consequently, 
the enhancement of the quark fluctuations is, in this model, mainly due to the amplification of the Boltzmann factor 
exp(— M/T) as the constituent quark mass is reduced. The suppression of Xqi Xi and Xuu above T c is due to the 
finite momentum cut-off of the NJL model. Indeed, the flavor-diagonal susceptibility Xuu in an ideal massless quark 
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FIG. 8: The diagonal Xuu susceptibility in an isospin symmetric system for /j, q = 0, 100, 200 and 305 MeV in the chiral limit 
normalized to Xuu C °' as a function of T/Tq. The calculations correspond to the vector coupling constants Gy = 0.3 Gs and 
G^/G^ = 1/3 and 1. 



gas with momentum cut-off A is given by 

V (free) /T 2 = ™ [* #P ^ T 
XUU ' ' W (2^( 1 + eP /T)2 

~ i _ l| e -A/r (1 + ^ 1 A2 

where the ellipsis in the second stands for terms that are exponentially suppressed for T < A. For A — > oo one obtains 
the ideal gas result Xuu^ /T 2 = 1, while for finite A the fluctuation are suppressed. At low temperature T« A, the 
correction terms in Eq. (|3.12| arc negligible, and x is independent of A. However, for temperatures on the order of A 
there is a strong dependence on the cut-off. With increasing temperature the suppression of Xqi Xi and Xuu grows 
stronger, as seen in Fig. [101 As seen in Fig. [10] this holds also for Xq an d Xi- However, the cutoff dependence of Xud 
is different. Here, since Xud oc Xq ~ Xli the cut-off dependence of Xq an d Xi partly cancel. Thus, the trend above 
and below To is different. At low temperatures the magnitude of Xud is reduced and at high temperatures enhanced 
with increasing A. We remind the reader that both parameter sets (a) and (b) yield the same transition temperature 
T = 177 MeV at fi q = (see Table 1). 

In Figs. [H] and [5] we also show the dependence of Xuu and Xud on the choice of vector couplings. The dependence of 
Xuu on Gy is fairly weak, as seen in Fig. [5] This is because in isospin symmetric matter (pLj = 0) Gy contributes 
only to x/, which is much smaller than Xq ( c f- Fig. [4]). Consequently, Xuu is essentially determined by Xqi which 
depends on Gs and Gy\ 

A comparison of Fig. [S] and shows that the off-diagonal susceptibility Xud is much smaller in magnitude than \uu- 
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FIG. 9: The off-diagonal x-ud susceptibility in an isospin symmetric system for [i q = 0, 100, 200 and 305 MeV in the chiral limit 
normalized to Xuu^ as a function of T/T Q . The calculations correspond to the vector coupling constants Gy = 0.3 Gs and 
G ( y VG^ S) = 1/3 and 1. 



Still, Xud is an interesting observable, which may be used to identify the transition point. This is particularly the case 
at finite /i q where the off-diagonal susceptibility is changing sign when crossing the critical temperature. This behavior- 
is consistent with the recent LGT findings which shows neg ative values of Xud below and above deconfincment for 
f/, q = 0. At finite [i q and at T = T c the LGT results show [26| an abrupt change of \ud from negative to positive value. 
At the tricritical point, Xud, being proportional to Xq ~ Xii diverges as \q (see Fig- HI (d)). Also seen in Fig. [5] is a 
rather strong variation of Xud with the strength of Gy /Gy ratio. Above the chiral phase transition, Xq and xi are 

equal for Gy^ = Gy . Hence, the fact that in LGT Xud is very small above To, may be interpreted as a signature of 
the universality of the isosclar and isovector current-current interaction in the chirally restored phase. For fi q = 0, Xud 

vanishes also in the chirally broken phase if Gy = Gy . The fact that in LGT Xud is negative in this temperature 

range, is consistent with Gy^ > G^y \ as expected for baryons from large N c arguments. 



3.3. Susceptibilities at finite current quark mass 



So far, we have computed the susceptibilities in the chiral limit, i.e., using an effective Lagrangian with an exact 
symmetry, the chiral symmetry. However, the chiral symmetry of the QCD Lagrangian is only approximate due to 
the finite current quark masses. In the following, we account for the explicit chiral symmetry breaking and explore 
its influence on the quark number fluctuations. This study could be also interesting from the perspective of recent 
LGT results [26| where the susceptibilities were calculated for finite and large quark masses. 

In Fig. [TT]we show results obtained with the NJL model for the flavor diagonal Xuu and off-diagonal Xud suscep- 
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tibilities for different values of the current quark mass. These calculations were done both for vanishing and finite 
chemical potential as well as for several sets of the vector coupling constant. 

A comparison of the results for different quark masses clearly shows that Xuu exhibits a peak structure at finite 
temperature, indicating the phase change, for all values of m. However, at finite m the Xuu is a smooth function 
everywhere, whereas in the chiral limit it exhibits a non-analytic structure at the phase transition. This behavior 
indicates that the second order transition at m = is converted into a smooth cross-over at finite m. The above is 
even more transparent at finite fi q where the discontinuity of Xuu at T c disappears at finite quark mass. The behavior 
of the NJL model at finite m is in accord with recent LGT finding in 2-flavor QCD, which also shows a cross-over 
transition at finite quark mass. We note, however, that the second-order transition expected in the NJL model in the 
chiral limit is still not confirmed by LGT calculations. Some recent results even suggest that in 2-flavor QCD this 
transition could be weakly first order [54| . 

With increasing quark mass tli6 peak position of \uu 

/T 2 is shifted towards larger T, both at vanishing and at finite 
fi q . The height of the peak also depends on m and decreases with increasing current quark mass. The shift of the peak 
position in Xuu indicates an increase of the pseudo-critical temperature with increasing quark mass. Such an effect 
was observed in lattice calculations already some time ago [45| . Also the suppression of quark number fluctuations 
with increasing quark mass is found in recent LGT calculations both at vanishing and at finite [i q . This suppression is 
related to an upward shift of the hadronic mass spectrum [3l| , which leads to a reduction of the number of thermally 
excited hadrons ~ exp(— Mh/T). The suppression of Xuu in the NJL model is of similar origin; the number of thermally 
excited quarks is cut back by the corresponding increase of the dynamical quark mass. The increase of Xuu with \i q 
seen in Fig. [TT]can be understood in terms of the corresponding increase of the thermal factors ~ exp(/i 9 /T). Finally, 
at large temperatures T > Tq the NJL model shows a very weak dependence of the susceptibilities on the current 
quark mass independently of the value of jj, q . In lattice calculations of QCD the m g -dependence of thermodynamic 
quantities was also found to be weak for m q /T < 1. 

The influence of a non-zero current quark mass on the off-diagonal susceptibility Xud for different values of fi q and 
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FIG. 11: The diagonal and off-diagonal susceptibilities normalized to T for different values of the current quark mass as a 
function of T/Tq. For the pseudo-critical temperature for m = 5 MeV we find To = 179 MeV. The results correspond to fj, q — 
and 100 MeV and the vector coupling constants Gy = 0.3 Gs with Gy /Gy = 1/3 and 1. 



vector couplings is illustrated in Figs. [TT] (c)-(f). For [i q = and G v /G v < 1, Xud is finite and negative for all 
T, and approaches zero for large temperatures. A similar behavior is observed in 2-flavor QCD on the lattice [26| . 

When Gy /Gy is increased towards unity, Xud approaches zero at all temperatures. At finite fx q , the temperature 
dependence of Xud changes qualitatively. Below T c it is negative, while above T c it is positive, in qualitative agreement 
with the results of LGT calculations [26| . The off-diagonal susceptibility Xud is finite near the pseudo-critical transition 
for all values of m. 
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When discussing the influence of the model parameters on quark number susceptibilities we have allowed for vari- 
ations of the parameters. However, we have not considered a possible T and jl dependence of G$,Gv and A. It was 
recently argued that such a dependence is important for a quantitative comparison of NJL model results with LGT 
findings [551 ] - However, so far systematic calculations of the variation of these parameters with temperature and chem- 
ical potential are lacking. A possible temperature dependence of these parameters was obtained phenomenologically 
by comparing some observables with lattice results. 



4. SUMMARY AND CONCLUSIONS 



We have discussed the properties of quark number fluctuations within the framework of the Nambu-Jona-Lasinio 
(NJL) model. The model was formulated at finite temperature and chemical potentials for baryon number and isospin. 
In a mean field approach, we have shown how the fluctuations of different quark flavors are changing across the phase 
boundary. Such a study is of interest from the perspective of both heavy ion phenomenology and lattice gauge theory. 
In the first case we have explored the non-monotonic structure of net quark, diagonal and off-diagonal susceptibilities 
along the phase transition line. We have also discussed the critical region around the tricritical point in the context 
of heavy ion phenomenology. 

The results on different quark susceptibilities at finite quark mass are in qualitative agreement with recent findings 
on the lattice. Our study may give some insight into how lattice results may change when one approaches the chiral 
limit at vanishing and at large baryon chemical potential. This expectation is supported by the fact that the NJL 
model exhibits the same critical properties as one expects for QCD. 

Obviously, the NJL model differs substantially from QCD. This model does not contain all relevant hadronic degrees 
of freedom, which in QCD contribute substantially to the quark number susceptibilities. Moreover, the phase-space 
of dynamical quarks is suppressed by the ultraviolet cut-off. Consequently, the perturbative regime of QCD at high 
temperature and density is not reproduced by this model. Nevertheless, features of the quark number susceptibilities 
that probe the restoration of chiral symmetry can be studied in some detail in NJL model calculations. In particular 
the change in the behavior of the quark fluctuations near the critical point can be interpreted as an effective change 
of the vector interaction associated with the chiral phase transition. 
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APPENDIX A: ANALYTICAL RESULTS FOR THE PHASE BOUNDARY 



For massless quarks, the gap equation (|2.13p at the chiral transition can be obtained in closed form. At the second 
order transition, the gap equation has a non-trivial solution at M = 0, i.e. 

d 3 p 1 



1 =4G S N C Y, 



f 



(2tt) 3 p 



1 - n [ f +) (p, T, fif) - n) ] (p, T, flf) 



is satisfied. After some rearrangement one finds for jl u = fi c i = p 
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where the quadratically divergent term on the left hand side is due to the vacuum loop and I(T, jl) to the thermal 
loops. The explicit form of the latter is given by 
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where ^[z] = Y^n=r z n /n 2 is Euler's dilogarithm [561 ]. 

We now use (|A.3[) to explore the dependence of the phase boundary on the cut off A. A general variation of (j A.2|) 
yields 

ata 71-2 dl c . dl crT1 dl _ ,. A . 

+ 4GIN C SG ° = 9A + df 5T + aj S » ■ (A - 4) 



At fx = 0, dl/dfj, = by symmetry and 

9/ 



dA e^ A + 1 



A. (A.5) 



Thus, the requirement that the critical temperature at vanishing net quark density remains fixed, as in the right panel 
of Fig. [2j implies the following relation between the cut off A and the coupling constant Gs 

Assuming that the transition is second order everywhere, the shift of the critical value of the chemical potential at 
T = is given by 

This relation, which remains approximately valid also when the transition is weakly first order, provides a quantitative 
interpretation of the shift of the phase boundary shown in the right panel of Fig. [2] 

Using the fact that for reasonable parameter choices, e^ ±A *~ A ' <C 1 along the phase boundary, useful approximative 
expressions may be obtained for I(T, fi). Retaining the first two terms in the expansions of the dilogarithms and the 
logarithm in (IA.3[) . we find 

j(T, n) = — + — — T(T + A) (>(^- A ) + e -0O*+A) 



6 2 

+ jT(T + 2A) ( e 2 ^-A) + e -2/3(p+A)) _ (A-8 ) 

Within this approximation, one finds the critical temperature at /i = with an accuracy of ~ 10~ 5 . For non-zero ji 
the approximation becomes increasingly better, since the A-dependent terms drop out for T — > 0. We note in passing 
that a solution of Eqn. (|A.2|) is possible only if the left hand side is positive, i.e., for GsA 2 > ir 2 /(2N C ) ~ 1.64. For 
smaller values of the scalar coupling constant, chiral symmetry is never broken. 



APPENDIX B: DERIVATIVES OF EFFECTIVE CONDENSATES 



In this appendix, we summarized derivatives of the dynamical masses Mf and the shifted chemical potentials flf. 
These results are obtained from the gap equations (|2.13p ~ (|2.15[) by taking derivatives with respect to [i q and [ii as 

dM f dM f dM f , ,dfif dfif, 

Q = Q A f + O A f + Q a f + Q B f ' 

1 = dfiu dM u g(s) dM u c(v) 8M d C ( S) _ dMd_ c (v) 

OfXqJ dfigj U dfigj " dfiqj d 8/J.qJ d 

dflu (S ) ^Mn (V) , dp, d (S ) _ ^Mrf (V) 

± i = ^ + g (S) _ d M « C ( V1 dMd_ c (S) dM d g(v) 

+ |^ S > - + + |^^ V) , (B.l) 

OfJ,q,I OHqJ dflqj a OflgJ 



19 



where / ^ /' S {u, d} and the functions Af, Bf, C| S,V ' > and D ( ^' v ' are defined by 



i(S,V) 



4 +) (l-nW)+4" J (l-4- J ) 



,(+) 



(-)/ 



5/ = 4G S N C 



d 3 p -AI f 



Cf- V) =4G^ V) 7V, 



D (s,v ) = 4G (s I v) iVc 



(2tt)3 rs/ 

d 3 p — M/ 
(2tt) 3 

d 3 p 1 



„<+>(! -„<+)] 



n\ ' 1 — n 



(-)l 



(2tt) 3 T 



Solving Eq. (|B. 1[) one gets, 



dM d = dM u 

dfJ, q ,l d[J, q j 

dp, u 



B u (l + 2D d V ' S) ) ± B rf (l + 2Z?i V ' S ))l /F , 



9Hq,I 

djld 



= ± 



1 - (A u +A d )(l + 2D { ^ S) ) T 2B d (C^^ T cf' s) ) + 2^ V ' S) J /F , 
1 - (A, + A d )(l + 2D^^) + 2B 11 (Ci v ' s > T C7' S) ) + 2^ V ' S) 1 /F , 



with the function F expressed by 

F = 1 - (A lt + 40 [l + £>< S >(1 + 2D d v) ) + £< V >(1 + 2£>f ) + £>f + D ( d Y) 



2 



cf > + cf + c 



(V) 



B u Df - fl d D( s > 



C (V) _ C (V) 



1 + 2D 



(V) 



D 



(S) 



1 + 2^ 



B u D™+B d D™ C&+C, 



(S) 



The Xg and X/ susceptibilities are obtained by substituting Eq. (|B.3[) into Eqs. (|3.2p and 



(B.2) 



(B.3) 



(B.4) 



APPENDIX C: CRITICAL EXPONENTS IN LANDAU THEORY 



In this appendix we compute the mean-field critical exponents of the TCP in Landau theory. Since the fourth 
order term b(T, in (|3.6p vanishes at the TCP, it is necessary to include a sixth-order term in the expansion of the 
thermodynamic potential 



u(T,n 9 ,M) ~ to (T,fi q ) + ^a(T, N )M 2 + ^b(T,fi q )M 4 + ^cM e , 



(C.l) 



where M is the dynamical quark mass. We assume that c > and neglect its dependence on the temperature and 
chemical potential. At the TCP the coefficients a and b both vanish. Close to the TCP we retain the leading terms, 

a{T, = A(T - T TC p) + B(ji - fx TC p) (C.2) 
b(T, n) = C(T - T TC p) + D(ji - fiTCp) ■ 

Along the second order phase boundary, the coefficient a vanishes, while at the first order transition a > and b < 0. 
The dynamical mass is determined by minimizing the thermodynamic potential. The solutions of the gap equation 

dui 



dM 



M{a + bM 2 + cM ) = 



(C.3) 
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are given by M = and M 2 = — j- ± -^Vb 2 — 4ac. The quark number susceptibility is given by 

5 (a4) 

As the TCP is approached, a, b and M — > 0. Thus, the singular part of the susceptibility is given by 

#"-\w?TZ- <a5) 

When the TCP is approached along the second order phase boundary, a = 0. Consequently x q mg = B 2 /2\b\ ~ 
I Ma — A*tcp| — ■ In order to compute the critical exponent along the first order transition, we must first determine the 
relation between a, b and c along this path. At the first order transition there are three degenerate minima of the 
thermodynamic potential, one at M = corresponding to the symmetric phase, and two at finite M corresponding 
to the two realizations of the broken phase. Thus, we arc looking for a solution to the two equations 

aM 2 + \bM A + \cM 6 = (C.6) 
a + bM 2 + cM A = , 

for non-vanishing M. Such a solution, M 2 = — 3&/4c, exists when the coefficients satisfy the relation 16ac = 3b 2 . 
This relation defines the location of the phase boundary in the (T, /i q ) plane. Furthermore, it implies that the 
first-order phase boundary is asymptotically parallel to the second order one. Using the relation in (|C.5p . we find 
X™ 9 = B 2 /\b\ ~ \fi g — ^tcp] 1 - Thus, the critical exponent is identical to that along the second order line, but 
the pre-factor is twice as large. If the TCP is approached from the broken phase, along any path which is not 
asymptotically tangential to the phase boundary at the TCP, b 2 -C 4ac for points close to the TCP. This implies that 
the critical exponent is different, namely x q mg = B 2 /A^/ac ~ \fx q — [J.tcp\~* ■ 

Finally, we explore the scaling behavior of the dynamical quark mass near the 0(4) critical line and at the tricritical 
point. As discussed in the text, this behavior determines the critical properties of the quark number susceptibility. 
The non-trivial solution of the gap equation implies that for constant chemical potential M 2 ~ a ~ |T C — T\ at the 0(4) 
critical line, where b ^ 0. On the other hand, at the critical end point, where also b — > 0, M 2 ~ y/a ~ \Ttcp — T\z. 
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